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ABSTRACT 



An approximate theory for cavitation in pump impellers is 
developed on the basis of a one-dimensional potential flow 
model. The theory incorporates the unsteady term occurring 
in Bernoulli’s equation and accounts for blockage effects due 
to vane thickness. A comparison of the theoretical results 
with experimental data from mixed flow impeller tests shows 
good agreement near the condition of zero incidence. On the 
basis of the theoretical results a mechanism of cavitation 
breakdown is proposed. 



NOMENCLATURE 



A = pipe area upstream of impeller, ft. 
o 

A^ = area of flow passage at the impeller eye without vanes, ft. 

= vane blockage coefficient, dimensionless 

C = ^ — = vane coefficient, dimensionless 

V 1-B 

V 

C = shroud coefficient, dimensionless 
s 

C = meridional velocity, ft/sec. 
m 

D = impeller diameter, ft. 

g = gravitational constant, ft/sec. ^ 

net positive suction head, ft. 

K = cavitation index, dimensionless 

N = impeller rotative speed, rpm 

2 

P = pressure, Ib/ft. 

2 

= vapor pressure, Ib/ft. 

q = velocity vector, ft/sec. 

3 

Q = flow rate, ft /sec. 

t = time, seconds 

u = impeller wheel speed, ft/sec. 

Z = number of impeller vanes 
a = relative flow angle, radians 

6 = vane angle, radians 

. . -1 
C = vorticity vector, sec. 

X = blade thickness, in. 

(p = C /u, flow coefficient, dimensionless 
m 

2 

^ = velocity potential, ft /sec. 

SUBSCRIPTS 

0 = reference station upstream of the impeller 

1 = reference station outside vane row at impeller inlet 
= reference station inside vane row at impeller inlet 
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A THEORY FOR CAVITATION IN IMPELLERS 



INTRODUCTION 

It has long been recognized that the flow in t urboraachinery 
is by nature unsteady when viewed from a fixed reference frame. 
Dean [ij pointed out the necessity of unsteadiness by showing 
that Bernoulli’s equation, applied between the inlet and out- 
let of a pump, denies the addition of energy to the fluid 
unless the unsteady terms are included. In other words, energy 
cannot be added to a steady, irrotational flow of a friction- 
less fluid; only by making the boundary conditions unsteady is 
it possible to add energy to the fluid. 

As a consequence of the inherent unsteady flow occurring in 
impellers when viewed from a fixed reference frame, it is 
necessary to account for the unsteady effects in any theoret- 
ical attempt to predict cavitation or minimum pressure coeffi- 
cient. Accordingly, the object of this paper is to develop a 
theory for cavitation in pump impellers which accounts for such 
effects. Bernoulli’s equation is applied to the flow in the 
inlet region of an impeller and the unsteady term occurring 
therein is retained. An analysis is then developed to approxi- 
mate this term so that the final result yields the cavitation 
index as a function of the flow coefficient. 

One of the most surprising results of this analysis is that 
the unsteady term in Bernoulli’s equation is not of secondary 
importance but has a major influence on the pressure in the 
inlet region at values of the flow coefficient on either side 
of the zero incidence condition. Moreover, on the basis of the 



theoretical results obtained, a mechanism for cavitation 
breakdown at positive and negative incidence conditions is 
propo sed • 

The underpressure in the eye of an impeller has been 
considered to be dependent on two factors, the effect of 
the shroud and the impeller vanes. To obtain the total 
maximum underpressure in the eye, the vane pressure drop 
is added to the pressure drop caused by the shroud. Reck- 
oning in this manner, Spannhake and Wislicenus developed an 
expression for the cavitation index in the form 



H = 
sv 



Ki 

2g 




2g 




2 -, 



where Kj and K 2 are shroud and vane coefficients, respectively, 

Q denotes the flow rate and D denotes the impeller diameter. 
Gongwer [ 2 ] was able to draw a mean line through his experi- 
mental data by properly selecting values for Ki and K 2 and, 
in terms of the notation of this paper, his resulting empirical 
expression for cavitation index at head breakoff versus flow 
coefficient at zero incidence is 

K -1.4 + .085 11+ — rr 

° \ % 

Although this empirical equation represents Gongwer’s experi- 
mental data for impellers having similar eye designs fairly well 
it gives little insight into the mechanism involved and does not 
reflect the effect of such features as vane thickness or number 
and non-zero incidence angle of the flow. 



An analytical method to calculate the complete flow field 
in an impeller is outlined in Ref. . This stream filament 

technique was programmed on a large digital computer and used 
to design a number of mixed-flow impellers described in Ref. 

[6]. A comparison of the test results of these impellers 
with the theory showed the theory to have the proper general 
trends and give reasonably accurate results in the portion of 
the impeller where channel flow exists. However, in the inlet 
region of the impeller the theory gave pressures which were 
consistently too high, making it of little value for the pre- 
diction of cavitation. 

Cavitation in impellers and inducers has generally been con- 
sidered to be described by separated cavity flow at the leading 
edge of the vanes, similar to that which is observed in steady 
flow past a single foil. Analytical methods to describe this 
kind of cavitation breakdown [l] , [s] , [a] have been based on 
the model of a free-streamline flow through a cascade of flat 
plate hydrofoils. While this model is no doubt a reasonable 
representation of cavitating flow at large positive angles of 
incidence, it is of limited value near zero and at negative 
angles of incidence (overflows). The present theoretical 
findings indicate that at negative incidence a favorable pres- 
sure gradient exists and there is no tendency for a separated 
region to exist. 

One of the first rather comprehensive contributions to the 
literature on pump cavitation was published by Gongwer [ 2 ] . 

A series of impellers of approximately the same eye design 
with widely differing outside diameters were tested with the 



result that the effect of the differing outside diameter had 
no effect on cavitation. Thus, Gongwer concluded that the 
cavitation performance of an impeller is only dependent on 
the eye design and the angle of incidence that the fluid makes 
with the impeller vanes. 

A simplified flow model advanced by Stanitz [5j to describe 
the effect of vane blockage at zero incidence was applied by 
Wood, Murphy and Farquhar [bj to their experimental results 
from mixed-flow pump tests. This model accounts for the 
reduced area caused by blockage effects of the vanes according 
to the one-dimensional continuity equation and relates the 
pressure at a station upstream of the impeller to the pressure 
at a station just inside the impeller according to the steady- 
flow form of Bernoulli’s equation. The theoretical approach 
presented herein is essentially an extension of Stanitz^s 
model to include the time derivative term in Bernoulli’s 
equation . 

The results of two outstanding experimental investigations 
of cavitation in mixed flow impellers are presented in the 
companion papers. Ref. [6j and [sj . The tests carried out 
under carefully controlled conditions showed the effect of 
vane incidence angle, thickness and number as well as speed 
effects. Since the geometrical characteristics of the series 
of impellers tested are well-recorded ^ the results were found 
to be very valuable for comparison with the theoretical result 
developed in this report. 



THEORY 



Two probable types of cavitation are hypothesized: (a) 

localized cavitation at the leading edge of the vanes, (b) 
gross cavitation behind the leading edge, which will be 
considered to be the true cavi ta t i on, giving rise to a 
significant reduction in performance. While the blunt lead- 
ing edge part of the vanes is no doubt the critical point in 
the impeller, it is believed that cavitation here would be 
only a local phenomenon and probably would not result in a 
reduction in performance. Thus, we will consider the point 
just behind the leading edge as the critical point for gross 
cavitation and attempt a theoretical prediction of the pres- 
sure or cavitation index at this point. 

In essence, the basic problem is to calculate the minimum 
pressure coefficient on the impeller under unsteady flow 
conditions. Assuming the flow to be incompressible and 
inviscid, Bernoulli’s equation applied between some point 
in the inlet conduit, far enough removed from the impeller 
eye so that the flow can be considered locally steady, and 
some point within the impeller, gives: 



In this equation q is the local velocity measured with respect 
to a fixed reference frame. The subscript m denotes ’’meridonal" 
and "0” corresponds to the upstream reference point. The 
velocity potential ^ is defined as =V(^>. Neglecting 

the elevation term and adding and subtracting the fluid vapor 



p -P 



2 




5 



pressure, , gives 



P -P P-P 

o V , 

r 2 + 1 - T 



pc 



mo 



pC 



mo 



mo 



2 

C 2 
mo 



( 1 ) 



The classical theory of cavitation states that cavitation 
occurs when the local minimum pressure just equals the vapor 
pressure of the liquid. It is well-known, however, that in 
practice this simple concept is not always valid. Departures 
of the cavitation index from values calculated on the basis 
of this hypothesis are caused by so-called scale effects 
which arise from changes in size, velocity and fluid properties 
Fortunately, however, experimental evidence indicates [9^ that 
in the absence of thermodynamic effects the scale effects 
decrease with size and velocity so that for the pro to type , thi s 
hypothesis yields valid results. Thus, proceeding on the basis 
that the size and speed of the impeller are sufficiently large 
so that the condition of cavitation is reached when P=P , 

V 



equation (1) may be written 

2 






mo 



2 a<i> 

C 2 3t 
mo 



(2a) 



where 



P -P 

K = 1 ^ 2 + 1 

o h qC ^ 
mo 



(2b) 



is defined as the cavitation index. The cavitation index in 
this formulation is exactly equal to one minus the minimum 
pressure coefficient. Unfortunately, this definition differs 
by a factor of 1.0 with the index, a, which is in common use 
for uniform flow past bodies such as head forms and hydrofoils. 
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It is, however, consistent with the definition generally used 
for turbomachinery. 

According to equation (2) it is evident that in order to 
calculate the cavitation index, it is necessary to evaluate 
both the absolute local velocity, q , and at the critical 

a t 

location on the impeller. As a first approximation, q may be 
taken to be the meridional velocity at the point in question, 
which will be dependent on the shroud and hub configuration, 
corrected for vane blockage. The effect of vane thickness 
can be approximated by accounting for the reduction in flow 
area presented by the vanes according to one-dimensional flow 
principles. 

The term 3^/3t can also be approximated once q is known 
provided the kind of unsteady flow involved is recognized. 

That is, although the flow is unsteady when viewed from a 
fixed reference frame, it is steady in a reference frame rota- 
ting with the impeller. As a consequence, we may imagine a 
velocity potential field fixed to and rotating with the impeller. 
The term d^/dtj which denotes the time derivative of ^ measured 
in the fixed reference frame, may be calculated by observing 
the change in ^ as the ^-field rotates past a point fixed in 
space during a small time interval. 

If the vane thickness was zero and the vanes were aligned 
with the flow, the surfaces of constant ^ would represent disks 

or surfaces axisymmetric with respect to the impeller axis. 
d $ 

In this case would vanish since an observer fixed in space 

o t 

would observe only one value of ^ as these disks are rotated. 
However, when the vanes have a finite thickness or are operated 
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at some angle of incidence, the flow is disturbed and the 



surfaces of constant ^ are distorted so that they are no 



longer axisymme tr ic . In this case 



8 $ 



3 t 



would no longer be 



zer o . 



Figure 1 shows a schematic of the unperturbed flow where 
the vane thickness is zero and the vanes are perfectly aligned 
with the flow. In this case velocity potential 

lines are aligned with the vane velocity vector as indicated. 
For the flow depicted here, 3<l>/3t = 0 because the lines corres- 
ponding to constant values of ^ move in the direction of their 
orientation. An observer fixed in space would observe the same 
value of ^ as time progresses and the velocity potential field 
rotates. For this special unperturbed condition the cavitation 

index as given by equation (2) reduces to 

2 



K =(C - /C ) 

o ml/ mo 



This relationship which represents the lower limit of K^reflectS 

only the effect of the shroud on cavitation index since the vane 

have no thickness and are aligned with the flow. 

Figure 2 illustrates the case of finite vane thickness, i.e., 

a small perturbation of the flow depicted in Figure 1. For this 

case it is evident that C ^>C , where C ^ may be related to C , 

m2 ml m2 ml 

approximately by the one-dimensional continuity relationship 

C , 



ni-L 

^m2 1-B^ 



(3) 



where the vane blockage coefficient, , represents the ratio 
of the vane frontal area to the area of the passage without 



vanes. For purposes of evaluating the projected vane 

frontal width is taken as x/sin 6 where t is the vane thick- 
ness and 6 is the local vane angle indicated in Figure 2. 

9 ^ 

Since C ^ is greater than C , and C = t — , the velocity 
m2° ml m9x ^ 

potential lines must be more closely spaced in the axial 
direction for finite vane thickness. Nevertheless, the 
basic kinematic boundary condition on the surface of the vane, 
given by 

= u*n =-u sinS (4) 



where n denotes the unit normal vector to the blade and u 
denotes the local vane velocity vector, must still be satisfied 
as it was for the case of zero vane thickness. That is, in 
either case the normal derivative of the velocity potential 
must satisfy the condition: 



9n 



T = 0 



9 <^ 
9 n 



T= f ini t e 



=-u sinS 



From this condition, then, the new 1 and 2 lines for the 
per t urba t ed, f low having the same values as those indicated in 
Figure 1, can be drawn. Thus, in Figure 2 we have represented 
by dashed lines the new 1 and 2 lines which are more closely 



spaced in the x-direction as required since ^jjj2 ^^ml still 

9 ^ 

having the same value of on the vane surface. That is, to 

satisfy this normal derivative requirement the new 4 lines are 

constructed so as to intersect points B and C in Figure 2, 

9 9 , 

To evaluate A<l> can be written as A ^ Ax where 

Ax* is the axial distance indicated in Figure 2. The corres- 

A c 

ponding At is At= where Ac is the distance indicated in 

u 



9 



Figure 2 and u is the vane velocity. Thus, 



bu t , 



as 



3<l> 

3t 



At 



3^ 

9x 



Ax ’ 

Ac 

u 



indicated in Figure 2, 

A X ^ 

^ = tan 0-0 
Ac 



9x 



Ax ' 
Ac 



(5) 

( 6 ) 



for small values of 0. Thus, from equation (5) and (6) can 

be written, 



at 






( 7 ) 



3$ 

where has been replaced by C ^2 



3<l> 



By definition of the velocity potential, yy is equal to the 
velocity of the fluid in a direction normal to the vane surface 



From Figure 2 it is evident that this can be written as 

3<l> 



3<l> 

3n 



3x 



(-AxM 



An 



9^ Ax’ 

37 a7 



Considering the geometry of the triangle (BCD) , indicated in 
Figure 2^ it is evident that 



a= Ax tang 



(9) 



and 

Ax-Ax’=a tan0 ~ a0 (10) 

Eliminating a between equations (9) and (10) and substituting 
this result into equation (8) gives 



yy = - yy (l-0tan6) cosS (H) 

3 $ 

Substituting this expression for t — into the kinematic boundary 

0 n 

condition given in equation (4) yields 



0 =co t B“ 



C 2 
ra 



where, again. 



3x 



has been replaced by Substituting this 



expression into equation (7) yields the following expression 
for in terms of the vane velocity, angle and meridional 
velocity : 



IT ■ “ yr) 

m2 



( 12 ) 



Finally, substituting the relationship given in equation (12) 
3<l> 



for 



3 t 



and using the approximation equation (2) 



yields the following relationship for th'e cavitation index: 

C 

E 

"mo / '^mo ^ml ^ '"m2 



K = 
o 



+ 2 ^ 



(cotB- ^ ) (13) 



Since, in general, the upstream reference velocity, C , 

mo 

and the local shroud velocity for zero vane thickness, C , , 

m i 

differ due to the effect of the shroud, it is convenient to 
introduce into equation (13) to obtain 



, \2 /c 



K = 



m 1 



■'m2 



o \C / \C , 
\ mo / \ ra 1 



+ 2 



^ m 2 

0 ml 






where (j) =C /u is the flow coefficient based on the upstream 
0 mo 

meridional velocity, C - /C can be called a shroud coeffi- 

m 1 mo 

cient since, from a one-dimensional flow point of view, it 
represents the ratio of the inlet pipe reference area to the 
local passage area if vanes are not present. ^jjj2^^ml 



called a vane coefficient since it represents the ratio of the 



local passage area without and with vanes, respectively. From 
one-dimensional continuity considerations, as given in 



equation (3 ) , 



m 2 X 

C = -z = — r =vane coefficient (15) 

V C , 1-B 

ml V 



and 



C A 

C = ^ = — 7-^ = shroud efficient 

s C A, 

mo 1 



(16) 



where B is the vane blockage coefficient and A /A_ is the area 
V o 1 

of the reference section in the inlet pipe to the local area of 
the passage excluding the effect of vane blockage. Using these 
definitions, the final expression for the cavitation index 
becomes 

2 2 12 1 

K =C C + 2 C cot3-^ -jr- (17) 

osv V b Vi/; 

0 s o 

Equation (17) relates the cavitation index to vane and shroud 

coefficients, the vane angle and the flow coefficient. The 

first term represents the dynamic depression associated with the 

velocity-squared term in the Bernoulli’s equation while the 

second two terms together represent the contribution from the 

unsteady term. The equation is valid at any radial location on 

the vane if S and (j>o are interpreted as the local vane angle and 

local flow coefficient, respectively. However, numerical evalua 

tion of the equation shows that the maximum value of K occurs 

o 

at the vane tip as might be expected. Accordingly, 3 and c|) 

0 

should be interpreted as the vane setting angle and ratio of 

to the vane wheel speed at the vane tip, respectively. It 
may be noted also that, although the ’’suction” side of the vane 
was used in the analysis, the resulting equation is the same 
regardless of which side is employed. It is only necessary that 



S be interpreted as the angle of the vane face in question 
if they are different. 

As indicated previously the cavitation index is, by 
definition, related to the minimum value of the pressure 
coefficient according to 

(min) 

where based on the pressure and velocity at the up- 

stream sub "o" reference station. A value of K =1.0 therefore 

o 

corresponds to the condition where the local pressure is equal 
to the upstream reference pressure. This, of course, would 
correspond to the minimum possible value of the cavitation 
index excluding effects of both the vanes and shroud. Values 
of less than unity would simply indicate that the local 

pressure was greater than the reference pressure and, these 
values are of no significants. 

The condition where the flow just inside the vanes at station 
’*2" is aligned with the vanes is of particular interest. This 
condition is defined by (j>2=t:an3 or, equivalently, = t an 6 / C . 
Substituting this value for <p^ into equation (17) yields simply 

K =C 2 (18) 

O S V 

because the second and third terms, corresponding to the contri- 
bution of 3$/8t, cancel each other. This result is essentially 
that of the blockage model proposed by Stanitz [ 5 ] , 

Values of (p which are less than tanB/C C correspond to 
underflows or positive incidence at the vanes and the contribu- 
tion from d(p/dt in equation (17) is negative, tending to decrease 
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. This effect, interpreted in terms of local pressure, 
corresponds to an increase in pressure and pressure coef- 
ficient . Although the measurements [6] show that the average 
pressure inside the vanes doeskin f act ^ increase for the 
positive incidence condition, it is also clear that the 
measured cavitation index increases, indicating a local 
depression. It is understandable then that Equation (17), 
which is based on one-dimensional flow considerations, 
disregards the local anomolies due to leading edge flow 
and predicts the increase in pressure at stations behind 
the leading edge. This positive incidence flow is dis- 
cussed in more detail and in relation to experimental 
findings in a later section. 

Flow at values of <j) which are greater than tanS/C C 

o s V 

corresponds to underflow or a negative incidence at the 
leading edge of the vanes. For this condition the last 
two terms in equation (17), which represent the contribu- 
tion from 3^/3t, give a positive contribution to . This 
effect, when interpreted in terms of pressure, corresponds 
to a depression and, consequently, reduced cavitation perfor- 
mance . 

COMPARISON WITH EXPERIMENTAL RESULTS 

Figures 3,4,5 and 6 show the experimental results from 
mixed-flow pump tests of Wood, Murphy and Farquhar [Sj and 
Wood [s] along with the present theoretical results obtained 
from equation (17). The experimental points represent values 
of the cavitation index as a function of the flow coefficient 
at the condition where the head rise curve just began to 



decrease due to cavitation. The details of the impeller 
geometries are given in Ref. and [^8j and, for purposes 

of application of equation (17), the features indicated in 
Table 1 were taken from these references. 







TABLE 1 


-IMPELLER INLET GEOMETRY 






Impel ler 




No . of 


Tan3 at Vane Tip 
























Vanes (suction (pressure 




B C 


C 








side) 


side) 


(Ave. ) 


V V 


s 


C 




4 


.291 


.321 


.306 


.0703 1.075 


1.04 


B 




5 


.299 


.329 


.314 


.088 1.095 


1.04 


A 




6 


.326 


.336 


.331 


.105 1.118 


1.04 


D 




8 






.46 


.1403 1.163 


1.04 


It 


may be noted 


that tanS in Table 1 


has two different values 


at the tip 


. This is 


due to the 


fact that 


the vanes were not 


o f 


constant 


thickness but increased 


in width 


wit h 


distance from 


the 


lead ing 


edge . 












On e 


part icular ly 


not ewor thy 


feature 


of the 


exp er imen t a 1 


data 


is the so- 


called velocity scale 


effect . 


The experimental cavita- 


tion index 


in Figure 


3-6 shows a 


decrease 


with 


impeller rotative 


speed with 


the higher 


speed runs 


showing 


the better agreement 


with 



the present theory. The authors of Ref.j^b] indicate that this 
velocity scale effect has an upper limit, however, since no 
distinct speed effect was observed in a comparison of the 4920 
rpm runs with some higher speed runs at 6500 rpm. it may be noted 
that this kind of dependence on speed has also been found in other 
experimental data for submerged bodies and, in general, the velocity 
scale effect has been found to decrease with velocity so that 
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eventually K approaches 1-C ^ . .In other words, the 

o p (min ; ’ 

cavitation index approaches the value as calculated on the 

hypothesis that cavitation occurs when the local pressure 

equals the vapor pressure. Since this classical hypothesis 

was the basis for the present theory the agreement with the 

high speed runs is encouraging. 

The three vertical line on Figure 3,4,5 and 6 are of 

special interest when comparing the present theory with the 

experimental results. The one line labled ({)^ = tan6 denotes 

the value of where the flow at the upstream reference 

station is such that the relative flow angle with the vanes 

would be zero. However, when the fluid encounters the vane 

leading edges it no longer has the velocity of C but has 

mo 

velocity due to the effect of the shroud. Thus, the 

condition of (f)j^ = tan6 (or (f)^= tan6/C^ ) has more significance than 

(p =tanS since the former indicates the condition where the 
o 

flow is approaching the vane leading edges at zero incidence. 
However, as pointed out by the authors of Ref, [b] , this condi- 
tion still is not the true condition of zero incidence, since 
due to vane blockage, the relative flow just inside the vane 
row would not be tangent to the vanes. The true zero-incidence 

condition occurs when 6« = tan6 or (p ==tan6/C C , For this condi- 

2 o s V 

tion equation (17) reduces to 




which correpsonds to the minimum value of K^, This value of 

(p is also indicated in the figures as a vertical line, 
o 



The authors of Ref.[6j have pointed out that their 
experimental data shows that the best cavitation perfor- 
mance, which corresponds to the minimum value of K , occurs 

o 

at the zero incidence condition of =tan3. Equation (17) 

corroborates this conclusion to the extent that K decreases 

o 

with decreasing values of (p . When (p^ = tan3 or (p =tanB/C C the 

2 2 ^ s V 

cavitation index becomes ^q”^s^v the second and third 

terms^ which arise on account of the 3<^/dt term in Bernoulli’s 

equation, cancel each other. However, for values of d> smaller 

o 

than tan6/C C , equation (17) shows a continued decrease in 
s V 

K while the experimental data show a minimum at <J) =tan6/C C 
o ^ o s V 

with a subsequent rise with decreasing values of . This 

continued decrease of K with decreasing values of <p simply 

indicates that the 3^/3t term occurring in Bernoulli’s equation 

gives rise to an increase in pressure over that existing at 

station 1. Although it is evident that could never drop 

2 

below (corresponding to cavitation outside the vanes at 

station 1) experiments show a rapid increase in below 
(p =tan3/C C . This rather extreme departure seems to 

0 S V 

indicate that some other mechanism is responsible for the 
cavitation breakdown at positive incidence angles. 

The difference between the theory and experimental data 
occurring at positive angles of incidence, that is, at values 
of (j) t an 3 / C ^ is not without an explanation. The one possible 
explanation which presents itself is that for positive incidence 
the pressure, as calculated according to the present theory, 
increases just inside the vane row. (This is born out by pressure 
measurements of Ref. [bj .) This pressure increase in the direction 
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of the flow (or adverse pressure gradient) would naturally 
tend to result in leading edge stall. However, for negative 
incidence the present theory shows a decrease in pressure as 
indicated by the rising value of with (j)^. This reduction 
in pressure inside the vane row would provide a favorable 
pressure gradient and, consequently, the result would be a 
well-behaved, attached flow. Thus, although at the outset 
there may seem to be no difference between positive and 
negative angles of incidence, as would be the case in steady 
cascade flow, the d^/dt term in Bernoulli’s equation gives 
rise to an adverse pressure gradient in the one instance and 
a favorable pressure gradient in the other. As a consequence, 
the resulting flow is apparently altogether different for the 
two possible incidence conditions. 
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SUMMARY AND CONCLUSIONS 



A potential flow theory has been developed to predict 
cavitation in pump impellers which takes into account the 
unsteady term occurring in Bernoulli’s equation. The theory 
has been applied to mixed-flow pump impellers and test results 
show good agreement near the zero incidence and overflow 
condition. As the relative flow angle of attack increases the 
theoretical and experimental results diverge due apparently 
to the predominance of the leading edge cavity in causing 
performance reduction. 

Figure 7 is an attempt to depict the salient features of 

cavitation breakdown in pump impellers. At (|) ^ = t ang / C^ C^ as 

indicated on the figure the unsteady term in Bernoulli’s 

equation causes no pressure r educ t ion, and equation (17) 

2 2 

yields simply K^=C^ C^ . This point is referred to as the 
point of zero incidence. For points to the right of this 
po int , c or r e spond ing to values of negative incidence, the one- 
dimensional theory, which disregards local effects at the 
vane leading edge, shows reasonable agreement with the 
experimental results. This agreement is apparently attri- 
butable to the flow stabilizing effect of the favorable 
pressure gradient caused by the 3<l>/3t term occurring in 
Bernoull’s equation. Points to the left of the point 
<^^ = tang/C^C^ correspond to values of positive incidence. 

The sharply decreasing versus curve with increasing 

incidence angle, a, indicates the existence of the adverse 
pressure gradients in the inlet region. This adverse pressure 
gradient apparently results in leading edge stall, a phenomenon 
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which would tend to block the flow passage and result in 
drastic performance deterioration. 

Although the present theory does not predict the minimum 
value of the cavitation index with acceptable accuracy it 
does shed some light on a very difficult and important problem. 
At overflow (or negative incidence) conditions the unsteady 
term in Bernoulli’s equation has been shown to provide a very 
important contribution to the pressure reduction in the inlet 
region and the theoretical results generally agree with the 
experimental data. At underflow the time derivative term 
contributes to an adverse pressure gradient. 




FIGURE I. VELOCITY POTENTIAL FIELD 

ROTATING WITH IMPELLER - 
ZERO VANE THICKNESS AND 
INCIDENCE ANGLE. 
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FIGURE 2. VELOCITY POTENTIAL FIELD 
ROTATING WITH IMPELLER. 
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FIGURE 3. CAVITATION INDEX VERSUS FLOW 
COEFFICIENT — IMPELLER A. 
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FIGURE 4. CAVITATION INDEX VERSUS FLOW 
COEFFICIENT- IMPELLER B. 
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FIGURE 5. CAVITATION INDEX VERSUS FLOW 
COEFFICIENT — IMPELLER C. 
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FIGURE 6. CAVITATION INDEX VERSUS FLOW 
COEFFICIENT- IMPELLER D. 
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FIGURE 7. CAVITATION INDEX VERSUS FLOW 
COEFFICIENT. 
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